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ABSTRACT 



Aims. Convective collapse, a theoretically predicted process that intensifies existing weak magnetic fields in the solar 
atmosphere, was first directly observed in a single event by Nagata et al. (2008) using the high resolution Solar Optical 
Telescope (SOT) of the Hinode satellite. Using the same space telescope, we performed a statistical analysis of convective 
collapse events. 

Methods. Our data sets consist of high resolution time series of polarimetric spectral scans of two iron lines formed in 
the lower photosphere and filter images in Mglb2 and CallH, spectral lines that are formed in the high photosphere 
and the lower chromosphere, respectively. We were thus able to study the implication of convective collapse events on 
the high photospheric and the chromospheric layers. The physical parameters from the full Stokes profiles were obtained 
with the MERLIN Milne-Eddington inversion code. 

Results. We found that in all cases, the event was accompanied by a continuum bright point and nearly always by 
a brightening in the CallH images. The magnesium dopplergram exhibits a strong downfiow in about three quarters 
of the events that took place within the field of view of the magnesium dopplergram. For each of the 49 events we 
determined the duration, maximum photospheric downfiow, field strength increase and size. We found event durations 
of about 10 minutes, magnetic element radii of about 0.43 " and 0.35 ", before and after the event, respectively, and 
field strengths of up to 1.65 kC 
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1. Introduction 

Time series of solar photospheric convection as seen in 
the continuum at visible wavelengths and the correspond- 
ing magnetograms show a highly dynamic quiet sun. 
Movements, accumulation, and separation of magnetic 
elements are clearly determined by the granular flows 
(|Berger fc Titld[l996l) . iStenflal (fl973li determined the field 
strengths of the coalesced magnetic concentrations by using 
the ratio between Stokes V profiles of two iron lines with 
roughly the same excitation potential and line strength but 
with different effective Lande factors. This line-ratio tech- 
nique revealed that the spatial-resolution-independent field 
strength is around 1-2 kG. This is signific antly more tha n 
the equipartition field strength of « 500 G (|Takeuchilll999i ). 
which is obtained by equating the known kinetic energy 
density of the co nvective flows with the magnetic energy 
density. This led iParkeil (|l978l ) to a possible explanation 
for the increase in fleld strength and the observed kG fields. 
This scenario is illustrated in cartoon form in Fig. [TJ 

Numerical simulations of the convective collapse pro- 
cess up to the destruction of the st able kG flux tube have 
been carried out in 2D and 3D (iHasaiil Il985t [Takeuchi 



1999 : [Grossmann-Doerth et al.l 19981: Sheminqva fc Gadun 



2000ft . The work by iGrossmann-Doerth et al.l (|19980 sug- 
gests that, in some cases, the downward traveling mass can 
"bounce of" the lower layers in the atmosphere and create 
an upward moving flow that could lead to the destruction of 



the compressed flux tube. Another des t ructio n mechanism 
was suggested by ISheminova fc GadunI (12000') and involved 
the cancellation of opposite polarity elements, destroying 
magnetic flux and permitting the disintegration of the flux 
concentration. 

After years of theoretical predictions and numerical sim- 
ulations, the flrst observati onal signs of the convectiv e col- 
lapse scenario were given bv lBellot Rubio et al.l (|200lh . who 
studied asymmetric Stokes V proflles of two photospheric, 
near-infrared lines at fixed spectrograph slit positions. They 
found redshifts in the Stokes V proflles, which correspond 
to the downfiow phase during a convective collapse pro- 
cess; however, shortly afterwards, they observed a decrease 
in the magnetic fleld strength, after reaching only 600 G. 
They also observed at the same time a strong, 3-lobed 
Stokes V proflle and argued that this is due to the appear- 
ance of an additional blueshifted Stokes V proflle. This can 
be explained by a sharp change from downfiow to upfiow 
along the line of sight and is in agreement with an upward- 
traveling shock front. Bellot Rubio et al. (2001) applied an 
inversion code based on the thin flux tube scenario and 
conflrmed that the proflles could be fltted very well with a 
discontinuity in the line of sight velocity as the result of an 
aborted convective collapse event. The first direct observa- 
tion of a convective collapse event by [Nagata et al. ( 200^ 
shows a magnetic network element undergoing a transition 
from a weak field strength of 400 G to a 2 kG flux tube. 
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Fig. 1 . Convective collapse scenario according to iParked 
(|1978f ). On the left the magnetic flux tube is sitting in the 
intergranular lane and is strong enough to reduce convec- 
tive heating, which will develop an adiabatic temperature 
gradient in its radiatively cooled atmosphere. The dotted 
lines represent the magnetic field and the horizontal line 
shows the r = 1 level and the Wilson depression inside the 
flux tube. The convective collapse process assumes an initial 
mass downflow, represented by the big dark arrow pointing 
downwards. The temperature of the downwards traveling 
mass increases adiabatically, whilst the surrounding atmo- 
sphere shows a much steeper temperature gradient. This 
leads to an enhancement of the downdraft. The flux tube is 
then evacuated and forced to compress as the gas outside 
the flux tube is pressing horizontally against it, represented 
by the thin arrows. The limit of the compression is reached 
when the magnetic field pressure balances the outer gas 
pressure. The right side shows the compression of the flux 
tube, which leads to the intensification of the magnetic field 
strength. 

The increase of magnetic field strength was preceded by a 
strong downflow of 6 kms~^. They were also able to link 
the strong-field magnetic element to the formation of a con- 
tinuum bright point. They concluded from these observa- 
tions that they were witnessing a convective collapse event. 
We present new observations from the Hinode Spectro- 
Polarimeter (SP) with the addition of chromospheric in- 
tensity in CallH and Mglb2 dopplergrams revealing the 
dynamics of the upper photosphere and perform a statis- 
tical analysis of 49 convective collapse events identified in 
several time series taken with the Hinode SOT. 

2. Data capture and processing 

We obtained 13 one-hour time sequences with the Hinode 
Spectro-Polarimeter (SP) each covering an area of 41 by 
4.8 arcsec with a cadence of 58 seconds and a pixel size of 
0.16 arcsec. The Fel 630.15 nm and Fel 630.25 nm spec- 
tral lines were used with a spectral resolution of 0.021 A. 
Simultaneous Narrowband Filter Imager (NFI) images in 
the wings of Mglb2 and CaHH filtergrams with the 
Broadband Filter Imager (BFI) were recorded. The Mglb2 
filtergrams were taken in both the blue and red wings at 
114 mA with an exposure time of 0.41 s and a pixel size 
of 0.16 arcsec. The CallH filtergrams were taken with an 
exposure time of 0.15 s and a pixel size of 0.11 arcsec. For 
more information on Hinode's Solar Optic al Telescope see 
iKosugi et all (|2007t ): [Suematsu et all ()2008[) : lShimizu et all 
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Fig. 2. Snapshot of a time sequence from the reduced data. 
From left to right, the coaligned images show the CallH 
intensity, the magnesium dopplergram, where black corre- 
sponds to downflows, the continuum intensity, the longitu- 
dinal and the transverse apparent magnetic flux (obtained 
by integrating over the circular and linear polarization pro- 
files and multiplying by a known calibration constant). The 
left side of the Mg I b2 dopplergram is cut off due to instru- 
ment misalignment and the restricted fields of view. 



(I2OOI : llchimoto et all (l2008h . We used the SolarSoft IDL 
routines sp.prep and fg.prep to process the data to level 1. 
The sp_prep procedure corrects the SP data for dark cur- 
rent, flat fleld, removes spectral and spatial drifts due to 
thermal instability, and merges the data from the two SP 
beams. The fg_prep routine corrects the flltergram (FG) 
data for dark current, flat fleld, and camera defects. 

We then spatially aligned the FG and SP data. Due to 
the small field of view in the narrowband filter and a fixed 
pointing offset between FG and SP data, part of the SP field 
of view is outside the narrowband filter field of view. The 
magnesium dopplergram therefore does not cover the whole 
SP field of view. The FG data was resampled to the SP tim- 
ing using a nearest-neighbor algorithm, with the central slit 
of the SP as a reference. The magnesium dopplergrams were 
obtained by dividing the difference between the blue wing 
and the red wing images by the sum. As the line-formation 
of the Mglb2 line is quite complicated and there is no 
simple response height that can be assigned to the mea- 
sured doppler signal (M.Carlsson, private communication) 
we could not obtain a calibration curve for the velocities 
in kms^^ from observations or simulations. However, the 
formation layer of the Mg I b2 dopplergrams was estimated 
by comparing calculated Mglb2 dopplergrams with snap- 
shots from MURaM simulations (|Vogler k. Schussleill2003l . 
see code description) , and the images revealed a strong sim- 
ilarity in the highest photospheric layers (N. Vitas, private 
communication) . In Fig. [2] we show one scan from a time 
sequence. 
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Table 1. Statistics of convective collapse events (see section 



Type of event 


max V down [kms ^] 


start B [G] 


peak B [G] 


r 1 [arcsec] 


r 2 [arcsec] 


duration [min] 


11 without Mg downflow 


3.80 


819 


1093 


0.48 


0.41 


10.36 


28 with Mg downflow 


3.92 


886 


1356 


0.48 


0.34 


10.37 


10 Mg out of FOV 


3.33 


900 


1262 


0.34 


0.31 


9.09 



We used the Milne- Eddington inversion code MERLIN, 
developed by the Community Spectro-polarimetric 
Analysis Center Q, to retrieve the ma gnetic field parame- 
ters from the observed Stokes profiles (jSkumanich fc Litei 
llQSTtlLites et al.ll2007D . 



with magnesium downflow 
no magnesium downflow 



3. Results 

Convective collapse events were identified by eye in the ob- 
served data and the inversion parameters. Each time series 
was checked for sudden substantial increases in magnetic 
field strength preceded by a strong photospheric downflow. 
We found that in all these cases, the event was accompa- 
nied by a continuum bright point and nearly always by a 
brightening in the CallH images. The magnesium dopp- 
lergram exhibits a strong downflow in about three quarters 
of the events that took place within the field of view of the 
magnesium dopplergram. 

Table [T] shows the statistics of the identified events in 
the 13 hours of repeated 4.8 x 41 arcsec^ raster scans. To 
obtain the properties of the magnetic elements we followed 
them by eye from one time step to the next. The max- 
imum of the line-of-sight velocity and the magnetic field 
strength are averages over the three highest data points, 
where the starting magnetic field was taken at the point of 
the strongest downflow. The majority of the events show 
equipartition field strengths at about 130 seconds before 
that, but it is difficult to trace the magnetic element pre- 
cisely at that time. The size was determined by tracing the 
magnetic element in the magnetograms by eye before (ra- 
dius 1) and after (radius 2) its collapse. The duration of the 
event was determined by taking the time when the maxi- 
mum downflow occurred as the start point until after the 
magnetic field had reached its maximum and declined to 
values comparable to equipartition field strength. 

The values in Table [1] are averages, which might be 
misleading were the distribution strongly asymmetric. The 
histograms of the line of sight velocity and magnetic field 
strength are therefore shown in Figures [3l [4] and \5\ These 
histograms indicate that the distribution of events with 
downflows seen in the high photosphere are similar to those 
obtained from the 11 events that do not exhibit downflows. 

Figure [6] shows a typical sequence in the calcium inten- 
sity, magnesium dopplergram, where black corresponds to 
downflows, the continuum intensity and the photospheric 
line of sight velocity. Clearly, a bright point can be seen 
in the continuum as well as in the Ca II H images forming 
just after a strong photospheric downflow. The magnesium 
dopplergram shows a strong downflow, too. The flux tube 
is therefore evacuated not only in the low photosphere but 
also in the high photosphere. As shown by the statistical 
analysis, this is the case for most events. An example of an 
event without significant downflows in the Mg I b2 doppler- 
gram is shown in Fig. [T] 
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Fig. 3. Histogram of the starting values of the magnetic 
field strength in the magnetic elements before the convec- 
tive collapse showing the number of occurrence divided by 
the total number of events. The dashed lines correspond to 
events without a downflow in the magnesium dopplergram. 
The solid line represents the 28 events exhibiting downflows 
in the magnesium dopplergram. 
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^ see http://www.csac.hao.ucar.edu 



Fig. 4. Histogram of the maximum field strength reached 
after the intensification in the same format as Fig. [31 



The evolution of a case with magnesium downflow sig- 
nals (c.f. Fig. [6]) is best seen in Fig. [H where we plot the 
photospheric line-of-sight velocity, magnetic field strength 
in the center of the magnetic element and the bright point 
intensity in the CallH filter. There seems to be an oscil- 
latory behavior in each parameter. This is seen in most 
observed events and will be the subject of future research. 
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Fig. 6. A convective collapse event asso- 
ciated with downflows in the upper pho- 
tosphere. The field of view is 4.8 " by 
4.8 ". The columns show the CallH in- 
tensity in a linear scale, the magnesium 
dopplergram in arbitrary units, the con- 
tinuum intensity in a linear scale, and the 
line-of-sight velocity derived through in- 
version of the Fe I line profiles, with posi- 
tive values corresponding to downflows. As 
discussed in section [5] we do not have a 
calibration curve for the magnesium dopp- 
lergram. The units are therefore arbitrary 
and black corresponds to downflows in the 
magnesium dopplergram. Green contours 
indicate the location of line-of-sight flux 
directed toward the observer, while red 
contours indicate line-of-sight flux directed 
away from the observer. The contour levels 
are at ±30 Mx/cm^ in the apparent longi- 
tudinal flux density. The blue contours in- 
dicate the transverse apparent flux density 
at 140Mx/cm^. 
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Fig. 7. A second example of a convec- 
tive collapse event, in the same format as 
Fig. 1^ This example does not show a sig- 
nature in the Mg I b2 dopplergram. 
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Fig. 5. Histogram of the maximum photospheric line of 
sight velocity in the same format as Fig. [3l 




500 E 



time [min] 

Fig. 8. Temporal development of the Ca H II intensity, the 
line of sight velocity and the magnetic field strength of the 
magnetic element from Fig. [5] Clearly, the maximum down- 
flow takes place shortly before the increase in the CallH 
intensity and the magnetic field strength. The event dura- 
tion in this case was longer than average, lasting about 15 
minutes. We show here the first 10 minutes, in which the 
oscillatory behavior in each parameter is clearly evident. 



4. Discussion 

As seen in Table [U the duration of convective collapse 
events is about 10 minutes. We can compare the duration of 
our events, i.e., the timespan from the moment of maximum 
downfiow until the time that the field has again declined to 
equipartition strength, to the lifetim e of magnetic el e ments , 
which has been determined by, e.g., Ide Wiin et al.l (|2008f) . 
They applied an automated feature-tracking algorithm to a 
six-hour sequence of magnetograms taken with the Hinode 
SOT, and found that magnetic elements on average have a 
lifetime of 10 minutes. They note that magnetic elements 
intermittently drop below the detection limit of the algo- 
rithm. It is tempting to speculate that the similarity be- 
tween their element lifetime and our event duration results 
from the ability to track elements only while they are in a 
collapsed state. 



Another process showing similar observational sig- 
natures is the emergence of magne tic flux. Th i s pro - 
cess has been well a n alyzed by, e.g., De Pontieu (I2002D: 
'Ce nteno et al.l (120071); [Martinez Gonzalez fc Bellot Rubiol 
f|2009D . ICenteno et al.l ^2007]) witness two opposite polar- 
ity patches emerging within a granular structure just after 
a horizontal field appears. We could be observing an emer- 
gent growing magnetic field instead of the enhancement of 
a field by the convective collapse process. However, the dif- 
ference is that the events we examine take place in the 
intergranular lanes, and we can observe in most cases the 
accumulation of fiux in the intergranular lane just before 
the process takes place. In the event shown in Fig. [SI we 
can clearly see the development over minutes before the 
event takes place. The positive polarity fiux (green con- 
tours) is carried by the granular fiows and has been present 
in the intergranular lane during 20 minutes before the sud- 
den increase of magnetic field. The opposite polarity patch 
near the observed event has been growing also through ac- 
cumulation and has been also present in the intergranular 
lane for several minutes. It is therefore unlikely that we are 
observing a fiux emergence event. The event in Fig. [7] does 
not have an opposite polarity near it before or during the 
collapse. 

It is still unclear what mechanism causes the de- 
struction of the st r ong m agnetic elements. According to 
iBellot Rubio et ahl (|2001h . an aborted convective collapse 
event can cause upward traveling shocks, which manifest 
themselves as asymmetries in the Stokes V profiles. The 
events observed here could fall into this category and ac- 
tually represent unfinished convective collapse events, or 
could be weakened or destroyed by upward flows after the 
collap se, as seen in simulations bv lGrossmann-Doerth et all 
()l998l ). Taking a closer look at the development of the 
Stokes V profiles (Fig. [9]) of the event shown in Fig. [6l 
we observe a redshift caused by the downfiowing material, 
but also strong asymmetries and 3 lobes in Stokes V pro- 
files in the center of the magnetic element. They could be 
signs for strong velocity gradients along the line of sight. 

Figure [10] shows two examples of the observed asym- 
metric Stokes V profiles, which can be either caused 
by several magnetic concentrations with different veloc- 
iti es in the same resolution element, or, as shown by 
.Bellot Rubio et al] ((1997), magnetic elements with strong 
velocity and magnetic field gradients along the line of sight. 
In the first case the net circular polarization (NCP) should 
be zero as the individual Stokes V profiles are antisymmet- 
ric. In the second case there should be a substantial NCP 
(more than 2 percent). We determined the NCP for several 
asymmetric Stokes V profiles, and find that all of them are 
around or below 1 percent. We therefore cannot confirm 
these events as being aborted convective collapse events or 
upward-traveling shocks. 

However, we did find examples consistent with another 
process that may weaken the magnetic elements. Figure ITT] 
shows magnetic flux of the opposite sign near to our in- 
vestigated collapsed magnetic element. The fiux is carried 
by convective motions toward the magnetic element. Then, 
horizontal flux is observed, and the inbound opposite polar- 
ity flux disappears almost completely. At the same time the 
field strength decreases sharply within the collapsed mag- 
netic element. This flux cancellation scenario occurs several 
times in our datasets and might be the leading cause for the 
destruction of the newly formed kG fields. This is in agree- 
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Fig. 9. Temporal evolution of the normalized Stokes V pro- 
files from the event in Fig.[6l The snapshots are 116 s apart 
and the time runs from left to right and top to bottom. 
The vertical lines correspond to the position of the inten- 
sity minimum averaged over the whole map. We chose the 
pixel at the center of the magnetic element. The strongly 
distorted Stokes V profiles are seen after the maximal down- 
flow and until the maximal field strength has been reached. 
This justifies the use of a Milne Eddington code to obtain 
the line of sight velocity and the starting and peak mag- 
netic field strength, as the Stokes profiles have a normal 
shape at these points. We show in red the fit obtained by 
the inversion code for the Stokes-V profiles at the beginning 
and end of the collapse. 



ment with simulations bv lSheminova fc GadunI (|2nOC)h . who 
state cancellation as the most common process of flux dis- 
persal. As mentioned previously, we have to consider that 
we might actually be observing an emerging flux event. 
We do observe in the event of Fig. [11] opposite polarity 
flux emerging and growing in the granular structure in the 
vicinity of the observed magnetic element, but at that point 
the convective collapse event has already taken place and a 
bright point is already visible in the continuum. 

The Stokes profiles were fitted using a variable filling 
factor. As we are in the weak field regime, the product of 
the filling factor and the magnetic field strength is very well 
determined by the inversion, but the filling factor and the 
magnetic field strength on their own are significantly less 
well determined. We therefore also computed the magnetic 
flux density per pixel inside the magnetic element, i.e. fill- 
ing factor multiplied by the magnetic field, and averaged 
these values. Comparing the magnetic flux density before 
and after the collapse we found an increase in the vast ma- 
jority of events, which is in agreement with the convective 
collapse scenario. The exceptions may be a consequence of 
the difficulty in determining the correct size of the magnetic 
elements by eye. 
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Fig. 10. Two examples of strongly distorted Stokes V pro- 
files observed within a kG magnetic element. The colors 
denote the blue or red lobe of the profile. We substract the 
blue area from the red area and divide by the sum to ob- 
tain the net circular polarization (NOP). All profiles show 
a NCP of less than 1 %. 
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Fig. 11. Flux cancellation event leading to the weakening 
of the collapsed magnetic element. The figure is in the same 
format as Fig.[6l but without the photosphcric line-of-sight 
velocity. 

In theory, the total fiux of the magnetic elements (on 
the order of 10^^ Mx) should not be altered by the convec- 
tive collapse as the increase in magnetic field strength is 
balanced by the decrease in area covered by the magnetic 
field. We do, however, find a systematic decrease in mag- 
netic flux, which may indicate a problem with determining 
the flux, the field strength, the filling factor, a real decrease 
in flux, or a combination of all of the above. 



5. Conclusion 

We found 49 convective collapse events in the quiet solar 
photosphere, all showing a brightening in the continuum 
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intensity and reaching peak magnetic field strengths of up 
to 1.65 kG, and with almost all events showing a brighten- 
ing in the Ca II H intensity. The evacuation of the flux tube 
is demonstrated by downflows in the lower photosphere, 
which are also observed in the upper photospheric layers. 
Our observations provide very clear evidence of the ubiq- 
uitous occurrence of the convective collapse process. Only 
about a quarter of these events do not show a significant 
downflow in the high photosphere. We also found that the 
duration of the events from the downflow until the peak- 
ing and subsequent decrease of magnetic field strength is 
about 10 minutes. In addition we confirm flux cancellation 
as a promising candidate for the weakening of the magnetic 
field strength in these str ong field elements as predict ed by 
simulations described bv lSheminova fc: GadunI (|200Cll ). 
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